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A BOUND ON GRASSMANNIAN CODES 


A. BARG* AND D. NOGIN f 


ABSTRACT. We give a new asymptotic upper bound on the size of a code in the Grassmannian space. The 
bound is better than the upper bounds known previously in the entire range of distances except very large 
values. 


1. Introduction. Let Gfc n (M) be the Grassmann manifold, i.e., the set of A:-planes passing through the 
origin in M n . Our focus is the packing problem in Cf-.n , he., the problem of estimating the number of planes 
whose pairwise distances are bounded below by some given value 5, for a suitably defined distance function 
d(p. q ). This problem has attracted attention in the recent years for several reasons. As a coding-theoretic 
(geometric) problem, it is a natural generalization of the coding problem for the projective space PW" 1 
and a closely related case of the sphere in M”, both having long history in coding theory 0. This problem 
arises also in engineering applications related to transmission of signals with multiple antennas in wireless 
environment 0. Finally, 0 introduced a construction of Grassmannian packings which is closely related 
to the construction of quantum stabilizer codes, another subject of interest in recent years. 

There are several possibilities to define a metric on G /, vn 0. We consider the so-called chordal metric 
(projection 2-norm in the terminology of 0), which can be defined in two equivalent ways. By a well- 
known fact 0, given two planes p, q £ 6'/,.. n one can define k principal angles between them. This is done 
recursively as follows: take unit vectors x\ G p. y\ G q with the maximum possible angular separation 
and denote this angle by 9 In step i = 2 ,,k, take the unit vectors x G p,ii 1 {x\, ..., Xi- 1 ) and 
y G q,yi _L (yi,... ,y*-i} with the maximum possible angle between them and denote this angle by 
In this way we obtain the set of principal angles 0 < 9k < • • • < 0\ < V 2 ; moreover, [x \...., X}-) and 
(yi,..., y/.) form orthonormal bases in p and q, respectively. 

Let sin 0 = (sin 0 \...., sin0*.). For a matrix (vector) P let \\P\\ = Pp denote its Euclidean 

2-norm. Define the chordal distance between p and q as follows: d(p,q) = || sin 0||. It turns out 0 that 
the Grassmannian space with the chordal metric affords an isometric embedding in a sphere S r of radius 
r = \fk(n — k)/n in R( n ~ 1 )( n + 2 V 2 Xo describe it, let A p be a “generator matrix” of p, i.e., a A: x n matrix 
whose rows form an orthonormal basis of p. Then the orthogonal projection from M n on p can be written as 
lip = ApAp. Define a map <f> : G}., n —> S r as d>(p) = I\ p — k /nl n (the plane is mapped to the traceless part 
of the projection on it). For any p, the norm of d>(p) equals 11I ; , — k /nl n \\ = r. The main result of |]4]j is that 
the mapping ( I> is an isometry in the sense that 

(i) d 2 ( P ,q) = V2||n p -n f/ || 2 . 

We call a collection of M points in Gk. n with pairwise distances at least 5 an (M, 6) code in the Grass¬ 
mannian space and call 5 the distance of the code. By 0 such a code gives rise to an (M, y/25) code 
C C S r , so any upper bound on the distance of C gives an estimate on the distance of 6y, vn . In particular, by 
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Figure 1. Bounds on R(6) for k = 3 (top to bottom): upper estimates 7 ?lp |2], Rr 
(Theorem 1); lower bound Rqy 0. 


Rankin’s bounds (8) for any (M, S) code Q, 

f k(n—k) 


M 


X <• ) n M—l 
0 ^ \ fc(n-fc) 


if M < n(n + l)/2, 
if M > n(n + l)/2. 


These bound are tight in the sense that there exist codes that meet them with equality El El- However, in the 
majority of cases, particularly, for codes of the large size, direct application of bounds on spherical codes to 
codes in Gk, n gives poor results because the image of Gk, n on S r forms a very sparse subset of it. 

We will be concerned with asymptotic bounds on M for a given value of the code distance 6. Let 
R = R(6) = limsup n _ >00 (l/n) In M be the largest possible rate of a sequence of codes with distance 6 in 
Gk, n - It is possible to compute the volume bounds on R analogous to the Gilbert-'Varshamov and Hamming 
bounds of coding theory 01- Namely, it is proved in m that for all 0 < <5 < Vk, 

(2) 2?gv(<5) <R<Ru(6), 

where 

Rqv(S) = —kln(6/Vk), 


Rh(S) = — fcln 


1 - 1 / 1 - 


2k 


The upper bound was subsequently improved in 0 relying on Delsarte’s linear programming method in the 
form developed in 0. The result of 0 is as follows: 

R < RlpV) := k[( 1 + s) ln(l + s) — sins], 

where s = ( fc / 2 )(^ — l)- This bound coincides with the result of 0 for k = 1 (codes in the projective 
space) and can be viewed as its generalization. For k = 1, Riji(d) < Rh(8) for all 0 < 5 < 1. However, 
for greater k the inequality 2 ?lp(5) < 1?h(< 5) holds only for 5 close to Vk and thus the linear programming 
bound provides a better estimate of R only for large values of the distance. For instance, for A: = 2,3 the 
crossing point is 6 = 0.74, 1.31 respectively. We note that Rup(Vk) = 0 showing that the lower bound 7?qv 
is tight for 5 = Vk. In this note we establish an improved upper bound stated in the following theorem. 

Theorem 1. 
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A BOUND ON GRASSMANNIAN CODES 


3 


Clearly Rr(5) < Rh(8) for all 6 G (0, \f~k- ] and R]\ ( \fk.) = 0. Moreover, Rr(5) < Rlp(5) for most 
values of 5 except for values in a small neighborhood of \fk. The intersection point 5* of the curves Rr and 
f?LP is given in the following table. 

k 2 3 4 5 10 

8* 1.37 1.717 1.992 2.231 3.161 

The behavior of the bounds for k = 3 is shown in the figure. 


2. Proof. The proof combines the isometric embedding of Gk, n in S r with an application of Blichfeldt’s 
density method similar to the arguments of Rankin ( 8 ). The intuition behind this method is as follows. 
Consider an ( M,8 ) code Q C G^.n- Denote by B ,5 = B$(x) a metric ball in G „ with center at x. Open 
balls of radius 5/2 centered at code points do not intersect, so no point of Gk, n can be contained in more 
than one such ball. The idea is to extend the radius <5/2 to some radius p so that while one point can belong 
to several balls, we can control the way the balls intersect and use some type of the volume argument to 
derive an upper bound on M. This idea, first suggested by Blichfeldt, can be viewed as a precursor to the 
well-known Elias bound of coding theory (see e.g., i TOl p.61]). 

Formally this idea is developed as follows. Under the mapping <f> : Gk, n —► S r an (M, 5) code Q is 
mapped to a spherical code A with minimum angular distance 2 a, where 5 = s/2r sin a. Let (3 be the angle 
given by sin / 3 = s/2 sin a and let p = s/2r sin P/ 2 . We compute 


(3) 



Let p G Q be fixed and let q G Gk. n be a point (plane) whose principal angles to p are given by 8 = 
(0i,... ,9k)- Let d = || sin 01| be the value of the distance between p and q. Consider the function on Gk, n 
defined by 


T P (q) 


{ 2cos/3 i „2 
r 2 sin 2 p \P 

0 


d?) if d < p 
if d > p. 


In other words, a p (q) can be viewed as a “density” defined on the metric ball B p C C/ lVn with center at a 
point p G Q and radius p. It depends only on the distance to the center (is spherically symmetric). 

Let us project the sphere S r radially on the unit sphere in S G id " -1 )!" -1 " 2 )/ 2 and denote the image of the 
code A by C. Applying <h followed by the projection to the ball B p with center at p transforms it into a cap 
on S with angular radius (3 and center at x = (\/r) ( \ > {p) on the surface of the sphere. The linear radius of 
the cap equals P = 2 sin P/ 2 . Letting q be a plane at distance d from p, we observe that the distance between 
x and z = d'iq) equals s = \j2d/r. The function r induces a function er on this cap defined with respect to 
x by 


0x(z) 


cos (3 
sin 2 f3 


(P 2 - s 2 ) 


for s < P and cr x (z) = 0 otherwise. A point z can belong to several caps with centers at points of the code 
C. The following lemma, whose proof is included for completeness, is due to 0. 


Lemma 2. For any point z G S, its total density satisfies 

TrO) < 1. 

x eC 


Proof: Let C C S be a code with distance 5 and let z G S be a point. Denote by x\,... x m G C the code 
points whose distance to z is at most P and let <7|..... d m be the values of these distances. We have 

^ ^ n m n 2 

-m(m - 1 )5 2 < - Y Y - x ik? = Y { m X P ~ ( Xki ) } 

i= 1 j,k= 1 i=l j k 
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n m 


E{ m E4<-(E^)} 


i =1 3 =1 


J=1 


= m 




— m — 


+X( m X4- (X- 

3 i =2 j j 


'31 


Without loss of generality let z= (1,0, 0,... , 0). Since 


dj = (1 - Xj i) 2 + x 2 2 H-h x’j n = 2(1 - Xji) 


we obtain the inequality 


1 


j 2 9 n 

|m(m-l)d 2 (Xf) ~X(X ; 

i j *=2 j 


'3 l 


which implies 
(4) 


d 2 J — 4m ^ d 2 + 2 m(m — 1)<5 2 < 0 . 


Let a z = J2?=i Vxj (z). We have 


Then 


a z = ^^-(m p2 -^d]) 
sin p — 


1 + cos (3 
2 cos P 


-a. 


yy d 2 = 4m sin 2 d /2 — cr 2 sin ft tan (3 = 4 sin 2 d /2 ^m — 

Using this in © we obtain 

16 sin 4 d/ 2 (m — 1 / 2(1 + sec/?)cr 2 ) 2 — 16 msin 2 d/ 2 (m — !/ 2 (l + sec/3)a 2 ) 

+ 2 m(m — l)d 2 < 0 . 

This inequality reduces to 4m(l — a z ) > a z tan 2 (3 which implies the claim of the lemma. 


Therefore also for any point q e Gk, n 

^2 t p(q) < !• 

P eC 

Let m(B p ) be the total mass of the ball computed with respect to the density r. From the last inequality we 
obtain 


(5) Mm(B p ) < Vol (Gk, n ) 

where Vol(Gfc )n ) is the total volume of the space. Let n{B p ) = m(B p ) / Vo\(Gk^ n ) be the normalized mass. 
We assume that k < n/2. The volume form on (the open part of) G'/ iV „ induces a distribution on the simplex 
of principal angles © = {( 0 i_, 0k) : V 2 > #t > ■ ■ ■ > dfc > 0 } given by 

k 

u>k, n = K{k,n) JJ(sindj) n ~ 2fc (sin 2 9i — sin 2 9j)dJd\... d6k, 

i=l 1 

where K(k, n) is a constant chosen from the normalization condition f G u>k, n = 1 (see, e.g., ( 6 |). Then 

p{B p ) = / r(||sin 0 ||)u; fejn . 

J 0:|| sin#||<p 

Asymptotic evaluation of an integral very similar to this one was performed in 0. We state the result in the 
following lemma whose proof is analogous to 0 . 
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Lemma 3. Let k be fixed and n — ► oo. Then 

/ n \ nk+o(n) 

M {Bp) = (tt) 

Substituting the last formula in © and taking logarithms we obtain 

— In M < —A; In -^= + o(l). 
n “ Vk 

Finally, using © and noting that r —> \/fc as n —> oo, we obtain the bound of the theorem. I 

The result of Theorem 1 can be also extended to the complex Grassmannian space similarly to an ex¬ 
tension to this case of the bounds © in ©. These estimates can also be extended to the quaternionic case 
according to the results of (ID- 
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